Is early thermalization really needed in A+A collisions? * 
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In this note we review our ideas, first published in year 2006, and 
corresponding results, including the new ones, which show that whereas 
the assumption of (partial) thermalization in rclativistic A+A collisions is 
really crucial to explain soft physics observables, the hypotheses of early 
thermalization at times less than 1 fm/c is not necessary. The reason for 
the later conclusion is that the initial transverse flow in thermal matter as 
well as its anisotropy, leading to asymmetry of the transverse momentum 
spectra, could be developed at pre-thermal, either partonic or classical field 
- Glasma, stage with even more efficiency than in the case of very early 
perfect hydrodynamics. Such radial and elliptic flows develop no matter 
whether a pressure already established. The general reason for them is an 
essential finiteness of the system in transverse direction. 

PACS numbers: 25.75.-q, 2^.10.Nz 

1. Introduction 

The problem of thermalization of the matter in ultra-relativistic A+A 
collisions is, certainly, one of the central in this field of physics. Even the 
answer to the basic question as for possible formation of new states of mat- 
ter, such as the quark-gluon plasma (QGP), in these experiments depends 
on a clarification of this problem. 

The experimental data support, in fact, an idea of thermalization. There 
are many dynamical models: hydrodynamic one [I], hybrid model - at first 
hydro evolution, then hadronic cascade [2 J , hydro-kinetic approach - coher- 
ent use of hydrodynamics and kinetics [3], that describe the soft physics 
observables in quite satisfactory way. All of them uses an assumption of 
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thermalization: either complete (local equilibrium), or partial (with viscos- 

ity). 

The discovery and theoretical description of elliptic flows at RHIC [4J 
play especially intriguing role in the problem of thermalization. It was 
cleared up [5] that the anisotropy of transverse momentum spectra, ex- 
pressed in terms of V2 coefficients, can be explained and basically described 
within hydrodynamic model for perfect fluid. Since the geometry of initial 
state in non-central A+A collision is ellipsoidal-like, where the biggest axis 
y is directed orthogonally to the reaction plane, then the collective velocities 
will be developed preferably along the smaller axis x - in reaction plane - 
since the gradient of pressure is larger in this direction. It leads to larger 
blue shift of quanta radiated from decaying fluid elements in x-direction 
at the freeze-out stage. This explains the anisotropy of the spectra, but 
to reach a quantitative agreement with the experimental data one needs 
to use very small initial time, r ~ 0.5 fm/c, to start the hydro-evolution. 
If it starts at later times, neither the collective velocities, nor their (and 
the spectra) anisotropy will be developed enough. These results, in fact, 
brought the two new ideas: first, that QGP, at least at the temperatures 
not much higher T c , is the strongly coupled system - sQGP, and so is an 
almost perfect fluid, and, second, that thermalization happens at very early 
times of collisions. 

A necessity for the thermal pressure to be formed as early as possible 
appears also in hydrodynamic description of the central collisions at RHIC. 
If one starts the hydrodynamic evolution in "conventional time" Tj =1 fm/c 
without transverse flow, the latter will not be developed enough to describe 
simultaneously the pion, kaon and proton spectra. 

The interesting observation related to the same topic was found in Ref. [6] 
concerning RHIC HBT puzzle. It was shown that to describe successfully 
pion, kaon and proton spectra as well as the interferometry radii Ri and 
explain the observed "puzzle" , R ou t/R s ide ~ 1> one needs i) initial transverse 
velocities developed by the time Tj =1 fm/c and ii) positive r^—t correlations 
between time t and absolute value of transverse radius rx on the surface of 
the particle emission during the system evolution. In the case of positive 
correlations the contribution to R 2 out -, which is proportional to the duration 
of the emission squared, can be compensated and, therefore, R ou t/Rside ratio 
can be reduced. In Ref. [7] the connection between these two conditions, i) 
and ii), is established. Namely, in hydrodynamic models the positive — t 
correlations at SPS and RHIC energies can appear only if fairly developed 
transverse flow already exist at the time Tj =1 fm/c. 

At the first glance, all those results for central collisions just conform the 
conclusion coming from non-cental ones: the thermalization and pressure 
are established at times essentially smaller than 1 fm/c and because of this 
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the radial/elliptic flows are already developed by this time. The crucial 
problem is, however, that even the most optimistic theoretical estimates 
give thermalization time 1 - 1.5 fm/c [8], typically, it is 2 - 3 fm/c. The 
discrepancy could be even more at LHC energies. 

The way to solve the problem was proposed in year 2006, Ref. [7j. 
It was shown that the transverse collective velocities can be developed at 
the very early pre-thermal stage with even higher efficiency than in the 
locally equilibrated fluid. This result was developed and exploited in Refs. 

In what follows we stick to the main line developed in these 

papers. 



2. Early developing of flow and its anisotropy 

A developing of the initial transverse velocities at the pre-thermal stage 
in ultrarelativistic A+A collisions is related to the complex problem of for- 
mation and evolution of chromo- electric and magnetic fields and their 
quanta - partons at very early stage of collision processes. In this note 
we analyze the developing of the transverse flows in the two opposite cases: 
evolution of non-interacting partons and fields and hydrodynamic expansion 
of strongly interacting system - perfect liquid. 



2.1. Non-relativistic analytical results. 

We start from the non-relativistic analytical toy- models. Let us put 
the initial momentum distribution of particles with mass m to be spher- 
ically symmetric Gaussian with the width corresponding to the thermal 
Boltzmann distribution with uniform temperature To, and without collec- 
tive flow: v(t = 0, r) = 0. Also the initial spatial distribution of particle 
density corresponds to the spherically symmetric Gaussian profile with ra- 
dius i?o- Let particles just to stream flreely. Then, according to [J3], the 
collective velocities, which can be defined at any time t according to Eckart: 

^7(t,x;p), (1) 
mr 

are 

where r is the radial coordinate and the homogeneity length Xhom is equal 
in this case to Ro. For the boost-invariant situation with initially elliptic 
transverse profile (corresponding Gaussian radii are R x = \hom,x and R y = 
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^hom,y) [14] the similar results for transverse velocities can be obtained at 
small t: 

Vi(t,r) ~ -g-*- (3) 

So, the collective velocities, which are developed in some direction i, are 
inversely proportional to the corresponding homogeneity length squared, 
^homi- The pressure plays no role in this collisionless process. The collec- 
tive flow in the expanding system appears because of a deficit of particles 
moving inward as compared to the particles moving outward. It happens 
due to finiteness of the system: one has less particles at the periphery than in 
the central part (the corresponding deficit is determined by the homogene- 
ity length) , and the flow develops during the evolution since the momentum 
of the particles that move outward cannot be compensated by the momen- 
tum of the particles which move inward, even if at the initial time such a 
compensation takes place. 

It is worth noting that, at least for fairly small time interval, the flows 
©J © for collisionless expansion coincide with flows, which are developed 
in hydrodynamically expanding perfect fluid (|14J) at the same initial con- 
ditions. 



2.2. Partonic free- streaming 

Let us start from the initial conditions inspired by the Color Glass Con- 
densate (CGC) approach. According to the results of Ref. [15] the time 
ro ~ 3/A s (where A s = g 2 p = Aira s p and p 2 is the variance of a Gaus- 
sian weight over the color charges of partons) is an appropriate scale con- 
trolling the formation of gluons with a physically well-defined energy. In 
Ref. [TTJ the results found in [15] for gluon momentum spectra in spa- 
tially transverse homogeneous case are transformed for the transversally 
finite systems and longitudinally smeared distribution instead of S(y — rf) 
(ri = llnp^ ~ 0,y = tin P0+PL ). Then the initial local boost-invariant 

\ ' 2 t xl ' a 2 Po—PL ' 

phase-space density / at the hypersurface (t 2 — x^) 1 / 2 = ro takes the form 

mi 

f( x ,P)\r = dN/d 3 xd 3 p\ TQ = ^^(r m T cosh(y - t?))^ 1 

(exp( V / ^^/T e// )-l)- 1 ^, (4) 

where m e ff = <22A S , T e jf = 03 A s ; a2 = 0.0358, 03 = 0.465. For finite 
systems we approximate the transverse profile of the gluon distribution by 
the ellipsoidal Gaussian p(b, ry) = po(b) exp (— r 2 /R 2 — r 2 /R 2 ) defined from 
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Fig. 1. Collective velocity developed at pre-thermal stage from proper time To = 0.3 
fm/c by supposed thermalization time 7j = 1 fm/c for scenarios of partonic free 
streaming and free expansion of classical field. The results are compared with the 
hydrodynamic evolution of perfect fluid with hard equation of state p = started 
at tq. Left : central (b=0), right : non-central (b=6.3 fm) collisions. 



the best fit of the participant number density in the collisions with the 
impact parameter b=(b,0). Constant a\{b) in Eq. is found then by the 
comparison of the total gluon number with results of Ref. [15] . 

Let us suppose that an actual thermalization happens about the time 
Ti = 1 fm/c > tq and partons just stream freely between tq and tj (see 
details for the r-evolution of the phase-space density in Refs. 0[9]), then 
the system is suddenly thermalized at Tj with the pressure P(e;Ti,x) cor- 
responding to the lattice QCD results. Then for tq = 0.3 fm/c, that 
corresponds to A s = 2 GeV [15], we find for the energy density profile: 
e(r ,6 = 0;r = 0) « 0.09A 4 /# 2 , Gaussian width R{r ,b = 0) » 5.33 fm 
(R x = R y = R). To describe pion spectra and interferometry radii one 
should use [II] A 4 s /g 2 « 6 GeV 4 and so, e(r = 0.3 fm/c, r = 0) = e = 67 
GeV/fm 3 . For non-central collisions with b = 6.3 fm: R x = 3.7 fm, R y = 4.6 
fm. 

The results for the collective flows are presented in Fig. 1 (b=0 and b = 
6.3 fm). We see qualitatively the same effect of the transverse flow and their 
anisotropy in non-cental collisions at pre-thermal stage as in the analytic 
non-relativistic examples. 
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2.3. Field expansion 

Let us describe now the flows at the early pre-thermal stage in Glasma 
approach [16] which deals with evolution/rearrangement of the classical glu- 
onic field. Let us again simplify the situation and consider extreme case of 
the boost-invariant evolution of the free field after moment tq ~ 3/A s [15] . 
Then we actually come to the Maxwell theory with 4-potential in space- 
time with the pseudo-cylindrical metric ds 2 = dr 2 — T 2 drj 2 — dx 2 — dy 2 . 

Imposing the gauge A T = used within the CGC concept, the dynamical 
variables of the model are <3? = A^, A4 (i = x,y). We will assume that these 
potentials and all the observables do not depend on space-time rapidity r\. 
It turns out that A x and A y are related by condition d{Ai = following 
from the gauge- fixing. 

The space-time evolution of the potentials is determined by the Maxwell 
equations which read 

d 2 Ai + -d T Ai - AtA = 0, (5) 
r 

dl<& - -d T $ - A T $ = 0, (6) 
r 

where Ay = didi. 

The components of the energy-momentum tensor, which are needed for 
finding collective velocity, are defined by the well-known formula: 

T Hv = _p^X F v x + glWF^F Xa , (7) 

where F Xa = d\A a - d a A\. 

If T^ v are known, the collective velocity is determined in accordance 
with the Landau-Lifshitz definition from equations: 

u 2 = l. (8) 



u x T Xa u° ' 

Now let us focus on the initial conditions for field equations. We apply 
the Glasma- like ones |16j : 

d T Ai\ T0 = 0, $| T0 = 0, (9) 



E 



^0, H z ~F xy \ To ^0. (10) 



These equations mean that, at the initial moment tq, the electric and mag- 
netic fields are longitudinal, while the transverse components vanish. So 
there is no field flow at tq. 
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With those initial conditions we get (at rj = 0) 

1 / d r „* 











To) 









Tu\r - <r T ) = -^- J +-( — J . (11, 

T tx \ T0 = 0. (12) 

Let us choose the initial energy density profile e(rx) = eoexp I —-^i 

to be the same as in the previous case of partonic system for b=0. 

To provide this we use, accounting for (jlip . the following transverse 
profile for potentials at the initial proper time tq: 



d rT V\ T0 = V^r T f{r T ), d T $\ T0 = VT^T f{r T ), (13) 



where /(rr) = y/2e(rj>) and a is a constant. Since the potentials VP, <3? are 
real, a bounds are < a < 1. The observable quantities T tt (r, ry, cp = 0, 77 = 
0), T tx (T, rx, (f = 0, 77 = 0), and t>(r, ry) are independent on parameter a. 

Further, we will treat these relations (together with d T ^f\ T0 = and 
<3?| T0 = 0) as the initial conditions for our problem. 

The solution for the field is expressed through the field amplitude at the 
initial moment tq: 

*o(fcr) = \fef(k T ), $o(M = VT^f(k T ), (14) 

here 

00 

f{kT) = Jfir T )jQ{k T r T )r T dr T . (15) 


Then the time-dependence of amplitudes takes the form 

§(r,fcr) = ^ (k T )^Y^[Ji(k T T )Y (k T T) - J (kTT)Yx(k T To)}, (16) 

<£(T,k T ) = $ (fe T )^-[J 1 (A: r r )y 1 (A ;T r) - ^(^^(^0)]. (17) 

To restore the spatial dependence, the Bessel-Fourier transform should 
be applied: 

00 

■$(r,r T ) = r T J ^>(T,k T )Ji(k T r T )dk T , (18) 


00 

$(r,r T ) = J $(T,k T )J (k T r T )dk T , (19) 
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where an integration over the angle has been already carried out. It allows 
one to calculate the evolution of the momentum-energy tensor (|7|) and then 
find the transverse flows according to Eq. ([8]). The results are presented at 
the Fig. 1. One can see that they are similar to ones obtained for partonic 
free streaming. 



2.4. Hydrodynamic evolution 

We compare the developing of the transverse collective velocities at pre- 
thermal stage by the time Tj= 1 fm/c with the hydrodynamic evolution of 
perfect fluid having the hard equation of state (EoS) p = |e. The initial 
conditions for hydro expansion: the time tq and the energy density profile 
are supposed to be the same as in previous cases (subsections 2 and 3). 
One can see from Fig. 1 that the hydrodynamic expansion is relatively 
less effective for producing transverse flows and their anisotropy than the 
non-interacting evolution. 



3. Conclusion 

We considered a possible developing of the transverse flows in differ- 
ent extreme cases: the evolution of non-relativistic free gas, partonic free 
streaming, classical field expansion and hydrodynamic expansion of the per- 
fect fluid with the hard equation of state. All the results are similar, though 
in relativistic situation the transverse flows, as well as their anisotropy, 
are developed at pre-thermal, either partonic or classical field, stage with 
even more efficiency than in the case of very early perfect hydrodynamics 
with hard EoS. So, one can conclude that the transverse flows and their 
anisotropy in non-central A+A collisions, in fact, do not depend on the 
nature of the systems formed and complicated kinetics of non-equilibrium 
processes. They are defined mostly by the homogeneity lengths in the ini- 
tial system. However, that is crucially important, the transformation of 
these radial and elliptic flows into observed high effective temperatures of 
proton and kaon spectra and anisotropy of pion spectra becomes possible 
only if (local) thermalization happens. Otherwise, e.g., at a free streaming, 
the initial partonic momentum distribution and its symmetry properties 
are preserved despite the flow and its asymmetry in non-central collisions 
are developed. So, the thermalization is certainly required for description 
of the hadronic spectra and i?2 coefficients but it does not matter whether 
thermalization is early or not. 



flowsl printed on January 25, 2009 



9 



Acknowledgments 

This work was supported in part by the Fundamental Research State 
Fund of Ukraine, Agreement No. F25/239-2008; the Bilateral award DLR 
(Germany) - MESU (Ukraine) for the UKR 06/008 Project, Agreement 
No. M/26-2008; and the Program "Fundamental Properties of Physical 
Systems under Extreme Conditions" of the Bureau of the Section of Physics 
and Astronomy of NAS of Ukraine. The research was carried out within 
the scope of the ERG (GDRE): Heavy ions at ultrarelativistic energies - 
a European Research Group comprising IN2P3/CNRS, Ecole des Mines 
de Nantes, Universite de Nantes, Warsaw University of Technology, JINR 
Dubna, ITEP Moscow, and Bogolyubov Institute for Theoretical Physics, 
NAS of Ukraine. 



REFERENCES 



[1] L.D. Landau, Izv. Akad. Nauk SSSR, Ser. Fiz. 17, 51 (1953) ;U. Heinz, J. 
Phys. G: Nucl. Part. Phys. 31, S717 (2005); T. Hirano, Nucl. Phys. A 774, 
531 (2006). 

[2] S.A. Bass, A. Dumitru, Phys. Rev C 61, 064909 (2000). 

[3] S.V. Akkelin, Y. Hama, Iu.A. Karpenko, and Yu.M. Sinyukov, Phys. Rev. C 
78, 034906 (2008). 

[4] K.H. Ackermann et al. (STAR Collaboration), Phys. Rev. Lett. 86,402 (2001). 

[5] P.F. Kolb, P. Huovinen, U. Heinz, H. Heiselberg, Phys.Lett. B 500, 232 (2001); 
U. Heinz, P. Kolb Nucl. Phys. A 702, 269 (2002). 

[6] M.S. Borysova, Yu.M. Sinyukov, S.V. Akkelin, B. Erazmus, and LA. 
Karpenko, Phys. Rev. C 73, 024903 (2006). 

[7] Yu.M. Sinyukov, Act.Phys.Polon. B 37, 3343 (2006). 

[8] Z. Xu, C. Greiner and H. Stocker J. Phys. G: Nucl. Part. Phys. 35 104016 
(2008). 

[9] M. Gyulassy, Iu.A. Karpenko, A.V. Nazarenko, Yu.M. Sinyukov Braz. Journ. 
of Phys., 37, 1031 (2007). 

[10] Yu.M. Sinyukov, S.V. Akkelin, Iu.A. Karpenko, Phys. Atomic Nuclei, 71 
(2008) 1619; arXiv: 0706.4066; N. Armesto et al, J. Phys. G: Nucl. Part. 
Phys. 35, 054001 (2008). 

[11] Yu.M. Sinyukov, Iu.A. Karpenko, A.V. Nazarenko, J. Phys. G: Nucl. Part. 
Phys. 35, 104071 (2008). 

[12] W. Broniowski,M. Chojnacki, W. Florkowski, A. Kisiel, arXiv:0801.436 
(2008). 

[13] J. Vredevoogd, S. Pratt. larXiv:0810.4325l f2008L 

[14] S.V. Akkelin, Yu.M. Sinyukov, Phys. Rev. C 70 , 064901 (2004). 



10 FLOWSl PRINTED ON JANUARY 25, 2009 

[15] A. Krasnitz, R. Venugopalan Phys. Rev. lett. 84 (2000) 4309; A. Krasnitz, Y. 
Nara, R. Venugopalan: Nucl. Phys. A 717, 268 (2003); A 727, 427 (2003). 

[16] T. Lappi, L. McLerran, Nucl. Phys. A 772, 200 (2006); 

F. Gelis and R. Venugopalan, Acta Phys. Polon. B 37, 3253 (2006). 



